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We report on magnetotransport measurements of multi-terminal suspended graphene devices. 
Fully developed integer quantum Hall states appear in magnetic fields as low as 2 T. At higher 
fields the formation of longitudinal resistance minima and transverse resistance plateaus are seen 
corresponding to fractional quantum Hall states, most strongly for v — 1/3. By measuring the 
temperature dependence of these resistance minima, the energy gap for the 1/3 fractional state in 
graphene is determined to be at ~20 K at 14 T. 
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In the low magnetic field regime, the integer quan- 
tum Hall (IQH) effect of graphene is marked by an 
anomalous half-integer quantum Hall conductivity a xy — 
g s (n + l/2)e 2 /h, where n is an integer and g s = 4 is the 
Landau level (LL) degeneracy resulting from the degen- 
erate spin and valley isospin degrees of freedom. This 
anomalous quantum Hall conductivity led to the observa- 
tion of the filling factor sequence v — ±2, ±6, ±10 [U 0- 
Subsequently, new broken-symmetry IQH states, corre- 
sponding to filling factors v = 0,±1,±4 have been re- 
solved in magnetic fields of B > 20 T, indicating the 
lifting of the fourfold degeneracy of the LLs [3j [4] . These 
filling factors have been suggested to be the result of vari- 
ous novel correlated states mediated by electron-electron 
(e-e) interactions [5]. 

In the strong quantum limit, e-e interactions in 2- 
dimensional electron gasses (2DEGs) can lead to the frac- 
tional quantum Hall (FQH) effect [6], many-body corre- 
lated states where the Hall conductance quantization ap- 
pears at fractional filling factors. In recent investigations 
of transport properties in two-terminal high-mobility sus- 
pended graphene devices [7, 8 , a quantized conductance 
corresponding to the is = 1/3 FQH state has been ob- 
served, suggesting the presence of strong e-e interactions 
in this system. However, due to the inherent mixing be- 
tween longitudinal and transverse resistivities in this two- 
terminal measurement [9 , quantitative characterization 
of the observed FQH states such as the FQH energy gap 
is only possible in an indirect way [10] . Although multi- 
terminal measurements on suspended graphene samples 
have been reported previously pTj H2] . the mechanical 
[T3] or thermal instability [14 of these samples has pre- 
cluded even the observation of a fully- quantized IQH ef- 
fect. 

Recently, the improvement of graphene mobility up 
to 8 m 2 /Vsec has been reported for substrate-supported 
graphene devices fabricated on a single-crystal hexago- 
nal boron nitride substrates [15] . Multi-terminal trans- 
port measurements performed on such devices in mag- 
netic fields up to 35 T reveal several FQH states whose 
filling factors are mostly integer multiples of 1/3. The 



energy gaps of these states have been measured for 
v > 1, exhibiting an unusual hierarchy among these FQH 
states [16]. However, the characterization of FQH states 
with v < 1, notably the 1/3 FQH state, could not be 
reliably conducted in these samples due to inhomoge- 
neous broadening near the charge neutrality point. As 
stronger e-e correlations are expected for this lower den- 
sity regime, further experiments on cleaner samples are 
desirable. 

In this letter, we report on the measurement of multi- 
terminal IQH and FQH effects in ultraclean suspended 
graphene samples in this low density regime. Filling fac- 
tors corresponding to fully developed IQH states, includ- 
ing the v = ±1 broken-symmetry states and the v = 1/3 
FQH state are observed. The energy gap of the 1/3 
FQH, measured by its temperature-dependent activation, 
is found to be much larger than the corresponding state 
found in the 2DEGs of high-quality GaAs heterostruc- 
tures. 

Our suspended graphene devices are fabricated using 
the method described in references [5J [11] : Mechanically 
exfoliated graphene samples are deposited on SiC^/Si 
substrates. Electrical contacts are made using electron- 
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FIG. 1: (color online). (a) Resistivity of a suspended graphene 
sample (SI) versus carrier density induced by back gate (re- 
sistivity is displayed on a log scale). Inset: Atomic force mi- 
croscope image of a typical suspended device, (b) Landau fan 
diagram R xx (V g ,B) at T=7 K measured in sample S3. The 
dark blue lines indicate minima in R xx . The dashed lines with 
integers indicate the corresponding filling factors. 
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FIG. 2: (color online). (a) R xx and R xy measured as a function 
of filling factor v at B =2 T and T = 1.7 K for sample SI. 

is plotted for v > 0.(b) and R xy measured as a 
function of magnetic field B at the hole density of n — 3.2 x 
10 11 cm -2 in samples SI (solid line) and S2 (dashed line). 
Numbers above horizontal lines indicate the corresponding 
filling factors. 



beam lithography followed by the thermal evaporation 
of Cr/Au electrodes. The Si02 under the flake is subse- 
quently removed via a chemical etch in buffered hydroflu- 
oric acid. A typical multi-terminal graphene device with 
a lateral size ~3 /im, suspended 150-200 nm above the 
Si02/Si substrate, is shown in the inset of Fig. 1(a). 

We studied a total of three 4-terminal and one 6- 
terminal devices where better quality results were ob- 
tained in three of the samples, labeled SI, S2 and S3. 
The longitudinal resistance R xx and transverse Hall re- 
sistance R xy were measured as a function of the ap- 
plied gate voltage V g on the doped Si back-gate, which 
tunes the carrier density n in the graphene. The car- 
rier density is determined using Hall measurements and 
n{V g ) = B/eR xy (V g , B), where B is the applied mag- 
netic field. The gate voltage is limited to less than 10 V 
to avoid electrostatic collapse of the suspended devices. 
The initial mobility of as-fabricated devices is typically 



less than 1.5 m 2 /Vs, comparable to samples on substrate. 
Sending a large current density (~0.5 mA//im) through 
a device heats up the graphene samples up to ~ 400° 
C [17], which typically removes many of adsorbed impu- 
rities, resulting in an extremely sharp peak in resistivity 
p as n changes from electrons to holes (Fig. 1(a)). In our 
study, the samples exhibit mobilities ranging from 8-15 
m 2 /Vs at the temperature T = 1.7 K, where most of our 
data was taken. In these ultraclean suspended samples, 
the Shubnikov de Haas (SdH) oscillations, resulting from 
the quantized cyclotron orbits are observable at relatively 
low magnetic fields. Fig. 1(b) displays a Landau fan dia- 
gram where R xx is plotted as a function of n and B in the 
low magnetic field regime (\B\ < 1 T). In this diagram the 
SdH oscillation minima (later developing to the quantum 
Hall minima) appear as strips fanning out from the origin 
point B = 0, n = with the slope dn/dB = ue/h. These 
strips survive in fields down to 0.1 T in our samples, in 
accordance with the mobility of ~10 m 2 /Vs calculated 
from conductivity measurements. 

As B increases, the observed SdH oscillations fully de- 
velop into the IQH effect. Fig. 2(a) shows R xx and R xy 
of device SI as a function of filling factor v = ne/hB 
(tuned by V g ) at a fixed magnetic field B =2 T. A series 
of quantum Hall states, i.e., zeroes in R xx and plateaus 
in R xy , quantized to values h/ve 2 with integer filling 
factors v = ±2, ±6, are observed within the gate bias 
window. More IQH statescan be observed using a field 
sweep at fixed gate voltage. Fig. 2(b) shows R xx and R xy 
measured as a function of magnetic field at a fixed hole 
density of 3.2 x 10 11 cm -2 for SI and S2. At least two 
well-defined plateaus with values (h/6e 2 ) and (h/2e 2 ) are 
observed, while the v = 1 broken-symmetry IQH state is 
being reached (h/e 2 ) at 14 T. We note that the devel- 
opment of this R xy Hall plateau is not complete, mea- 
suring only ~95 % of the full quantization value. This 
deviation of the Hall resistivity from the expected quan- 
tization value may be attributable to the presence of 
non-ideal disordered contacts which can introduce a non- 
equilibrium population of edge states that perturbs the 
quantization for the small samples used here [18] . It has 
also been suggested that the proximity between the cur- 
rent leads and the voltage probes could short-circuit the 
hall voltage in such small samples [13] . 

As we move to the low-density regime corresponding to 
v < 1, the FQH effect starts to be detected for B >10 T 
in all three samples (SI, S2 and S3). Fig. 3(a) displays 
R xx and R xy plotted as a function of B at a fixed den- 
sity of 1.9xlO n cm -2 (dotted line) and 0.96xlO n cm -2 
(solid line) respectively. At higher density, we notice that 
R xy increases further and R xx shows an additional dip 
at the field corresponding to v « 2/3. Upon further de- 
creasing the density, lower filling fractions come into the 
observable window set by the maximum probing mag- 
netic field (14 T), and R xx develops even deeper local 
minimum at the field corresponding to v ~ 1/3 with a 
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plateau-like feature in R xy . However, similar to the v = 1 
broken- symmetry IQH state, the corresponding features 
in R xy are not fully quantized in this low density and 
high Hall voltage regime. In sample S3 (Fig. 3(b)), in 
addition to minimum around v « 1/3, another minimum 
is visible around v ~ 1/2, but there is no feature close to 
a 2/3 filling fraction. 

We note that, when scaled by the filling factor, these 
local minima of R xx are robust features at different mag- 
netic fields and densities. Fig. 3(c) and Fig. 3(d) show 
the R xx of SI and S2 as a function of filling fraction v at 
different B and n. The R xx traces exhibit local minima 
corresponding to filling fractions v ~ 1/3, 1/2, 2/3 and 
1. For all three samples, we observe strong minimum for 
v = 1 and 1/3. But v ~ 1/2 only shows up clearly for two 
samples S2 and S3 while v ~ 2/3 emerges only in SI, indi- 
cating that these features are rather fragile and are more 
sample dependent compared to v =1/3. For example, in 
Fig. 3(d), in addition to the IQH state at v = — 1 (fea- 
ture C) two other notable features emerge as relatively 
deep minima in R xx (A and B), located at filling factors 
v = —0.34 and v = —0.64 respectively. Further confir- 
mation on the nature of these states can be provided by 
means of a Landau fan diagram where R xx is plotted as 
a function of V g and B (Fig. 3(e)). From the slopes of 
the lines in marked A and B we estimate that features A 
and B follow v = 0.31 ± 0.02 and v = 0.66 ± 0.02 lines, 
respectively. We thus assign features A and B to be the 
minima corresponding to the 1/3 and 2/3 FQH states. 

The strongly developed minima in R xx for the 1/3 state 
now allows us to probe the energy gap associated with it. 
To quantify the energy of this FQH state, we measure the 
temperature dependence of R xx . Fig. 4(a) and (c) display 
RxxiYg) measured at a sequence of different temperatures 
at 14 T for SI and S2, respectively. The minimum corre- 
sponding to the 1/3 FQH state, it^ m , increases as T in- 
creases. An Arrhenius plot for the R™ 71 (Fig. 4(b)) shows 
an activation behavior indicated by \og(R x n7n ) ~ T _1 in 
the temperature range between 9-22 K for SI. At lower 
temperatures (2-8 K) \og(RlV j7n ) deviates from a simple 
activation behavior , turning into the slower temperature 
dependence expected for variable-range hopping behav- 
ior. Similar trends were observed in early experiments 
on GaAs 2DEGs T9]. 

From the high-temperature activation behavior, the 
transport gap AE of the 1/3 FQH state can be obtained 
using i?™ m ~ exp(—AE/2kBT) at a fixed magnetic field. 
Fig. 4(d) shows AE as a function of B determined from 
the line fits in Fig. 4(b) for sample SI. A similarly de- 
termined activation energy gap for the 1/3 FQH state of 
sample S2 is also included in this figure (For this sam- 
ple the R™ 71 shows an activation behavior in the tem- 
perature range between 4-13 K). For both samples, we 
obtained similar magnitude of the energy gap, increasing 
with increasing 5, as is expected. 

Naively, the magnetic field dependence of energy gap 
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FIG. 3: (color online), (a) R xx and R xy measured in sample 
SI as a function of B at fixed density of 1.9xlO n cm -2 (dot- 
ted line) and 0.96xlO n cm -2 (solid line), (b)^* and R xy 
measured in sample S3 as a function of magnetic field B at 
fixed density of 0.8xlO n cm -2 . Numbers with vertical lines 
indicate the corresponding filling factors. The straight hori- 
zontal lines indicates the expected R xy value for the 1/3 FQH 
state, (c) R xx as a function of filling factor at B = 12 — 14 T 
measured in sample S2.(d) R xx as a function of filling factor 
dX B — 12 — 14T measured in sample SI. The letters A, B 
and C with vertical lines indicate minima corresponding to 
filling factors 1, 2/3 and 1/3 respectively, (e) Landau fan di- 
agram R XX (V 9 , B) at T=1.7 K measured in sample SI. The 
dark blue lines indicate minima in R xx - 



of 1/3 FQH state would be expected to increase with 
~ \[B, considering the e-e interaction energy scales with 
the Coulomb energy scale e 2 /e/#, where e is the dielec- 
tric constant and Ib is the magnetic length proportional 
to 1/y/B. As an alternative scenario, however, the acti- 
vation energy gap could be linear in £?, if the nature of 
charged quasiparticle excitations of this FQH state are 
associated with spin-flips in skyrmion-like excitations of 
a spin-polarized FQH ground state [20]. Given the field 
range and error bars of the data points, we cannot rule 
out either scenario at present. We remark the linear fit 
(dashed lines) yields a negative y-intercept ~ 20 K, which 
can be interpreted as the broadening T of the fractional 
state. This value is slightly larger than fi/r ~ 4 K, the ex- 
pected LL broadening estimated from the scattering time 
r obtained from the mobility of the sample. However, the 
y/B-&t (dotted lines) yields a negative y-intercept ~ 60 K 
too large to be considered a reasonable T. 

Considering the LL broadening T ~20 K obtained from 
the linear fit, the intrinsic energy gap for v =1/3 is es- 
timated to be Ai/3 = AE + T « 40-45 K at B =14 T. 
It has been predicted that the 1/3 state is both spin 
and valley-isospin polarized in the SU(4) configuration 
space [21], and calculations for this fully polarized 1/3 
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FIG. 4: (color online) (a) R xx measured at different temper- 
atures at 14 T for sample SI. (b) Arrhenius plots of R™ 71 
versus 1/T at different fields at v — 1/3 FQH state for sam- 
ple SI. The straight lines are linear fits to the data at the 
high-temperature activation range, (c) R xx measured at dif- 
ferent temperatures at 14 T for sample S2. (d) The activation 
gap as a function of field for two samples: SI and S2. The 
dashed lines are fits to B and the dotted lines are fits to \[B. 

state have given a gap value of C e 2 /eZ#, with the nu- 
merical constant C ~0. 05-0.1 [21j[22]. For graphene in 
vacuum, the dielectric constant from in-plane dynamic 
screening is estimated to be e =5.24 [23]. The predicted 
gap is then in the range Ai/ 3 =26-50 K at B =14 T, in 
reasonable agreement with the gap measured in our ex- 
periment. For comparison, the observed disorder reduced 
Ax/3 m graphene is at least 3 times larger than that of 
the 2DEGs in the best quality GaAs hetero junctions in 
a similar field range [20]. We further remark that A x / 3 
obtained in this experiment is much larger than the gaps 
obtained for FQH states associated to v > 1 considering 
those gaps are obtained at B =35 T [16]. This compari- 
son thus suggests an unusual robustness of the 1/3 state 
in graphene, inviting further investigation to elucidate its 
microscopic nature [24] . 



In conclusion, we have measured the quantum Hall ef- 
fect in multi-terminal suspended graphene devices. Both 
the broken symmetry IQH and 1/3 FQH states emerge, 
as manifested by minima in R xx and nearly quantized 
plateaus in R xy . From activation behavior of the R xx 
minima, the energy gap associated with the 1/3 FQH 
state is measured. 
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